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local extrema. We suppose that these meta-stabilities may influence 
the order of the phase transition in compact QED. 
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1 Introduction 

In compact QED (CQED) there are two distinct phases. The weak couphng 
regime corresponds to the Coulomb phase with massless photons and has been 
investigated in very detail in the extremely successful program of perturbation 
theory. However, there is also a second regime realized in the strong coupling 
region, which has become an important research topic since the 1980s. Here, 
magnetic monopoles cause the creation of electric flux tubes leading to a linear 
potential and therefore to confinement. Although the confinement mechanism 
in QED seems to be rather different from the one in QCD, the investigation of 
the former gives some inference to the confinement in QCD, which is far more 
complex to investigate due to its non-abelian structure. 

In the weak coupling phase the monopoles are supposed to disappear but they 
can survive in meta-stable states. There has been intensive discussion whether the 
phase transition between these two phases is of first or second order. For a long 
time it seemed that this transition is of first order mi2j. Some later investigations 
suggested a second order phase transition [3 HI EHH]. Recent studies [3 |H] strongly 
suggest a first order transition again. 

The investigation of the order of the phase transition is of vital importance 
since the continuum limit can only be obtained in the case of a second order 
transition. 

2 QED on a Space— Time Lattice 

Here we give a short summary of compact QED on the lattice. The Wilson 
action ^9] of QED is given by 

Sw = /? 5^(1 -cos (1) 

X,fl<U 

A common choice for the definition of the field-strength is 

ea^F^, = e^, , (2) 

where the physical angle 6^j,y G (— 7r,7r] is obtained by splitting off n^u{x) G 
{—2, —1, 0, 1, 2} the number of Dirac strings penetrating a plaquette 

Of,u{x) = 9^i,{x) + 27mf,^{x) . (3) 

6^y{x) and n^^{x) are gauge-dependent contrary to 9^,^. In our simulations we 
use the definition due to Lang et al. (TUl 

ea^Ffj,^ = sin6f,u , (4) 

where the field strength is a continuous function of ^^i, and takes into account 
27r-periodicity. This definition is achieved by a variation of the Wilson action [TT] 
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and is therefore in agreement with the GauB law on the lattice. As will be 
described further below, definition (jH) allows to investigate continuous shifts of 
monopoles by arbitrary non-integer distances, a property which we need for our 
investigations. 



3 On Magnetic Charges 

As mentioned in the introduction an interesting feature of U{1) gauge theory is 
the appearance of magnetic monopoles. So we will shortly recall the fundamental 
relations for Dirac monopoles. 

Dirac considered a 3-dimensional ball surrounding a magnetic point charge 
yielding the following relation 

g = I pmd^x = VBd^x = i Bdf= i Ads. (5) 

J ball J ball J sphere J 9sphere 

The boundary of a boundary like ^sphere is zero, which implies that the vector 
potential diverges at a single point on the sphere for non- vanishing g. One is 
lead to the picture of monopoles connected by arbitrary lines with a diverging 
vector potential along them. These lines are the famous Dirac strings. Their 
position is gauge dependent. So, there are two kinds of divergences. First the 
"true" physical divergences of monopoles and the unphysical "gauge" divergences 
of Dirac strings. 

Furthermore, Dirac showed that the electric charge e is quantized by the exis- 
tence of a single magnetic charge g since the wave function of an electric charge 
transported along C is modified by 

exp / Adi] (6) 



c ' 

Transporting e around the boundary ^sphere of Eq. (0), the non-ambiguity of i\} 
demands 

I / Ads = ^ = 27rn. (7) 

^Osphere 

Setting h,= \ yields for the smallest (n = 1) charge quantum 

e^-. (8) 

4 Identification of Monopoles 

Dirac monopoles in compact U(l) theory can be identified according to DeGrand 
and Toussaint [ISI as follows. 
The plaquette angle 

Q^v{x) = e^{x) + e^{x + /i) - d^,{x + e^{x) e {-An, An] (9) 
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has no direct physical meaning since the action has a periodicity of 27r only. 
Furthermore, the sum of six plaquette angles enclosing a cube vanishes since 
each link of the cube is counted twice but in opposite direction. In the language 
of forms the cube is the exterior differential of the six plaquettes (respectively 
angles), which are the exterior differentials of their links. Due to the Bianchi 
identity ("the boundary of a boundary is zero") the sum of the six plaquette 
angles of a cube vanishes. 

Dirac strings can be closed or start and end at magnetic monopoles. Thus, the 
magnetic charge in a cube, respectively the number of monopoles in it is given by 
the sum of in-going minus out-going strings. It is an integer multiple of (7 = 27r/e. 
The corresponding four-current is given by 

imS^(^) = '^e^,,^A,np„{x), (10) 

where we used the (forward) lattice derivative A^f{x) = f{x + z/) — f{x). Ob- 
viously, this current is source-free (A^jm,/^ = 0). Identifying monopoles this way 
allows to count the number of monopoles in each cube only. For the identification 
of monopole positions inside cubes we have to choose another method. 

This second method, which is due to Ref. ^1] detects monopoles by measuring 
the magnetic fluxes leaving cubes. The distribution of the flux over six plaque- 
ttes surrounding a cube reflects the position of the monopole inside. Then, the 
monopole current is given by the dual Maxwell equations, which read in covariant 
notation 

9^F^i, jm,u (-^-^) 

d^F^. = (12) 

with the dual field strength-tensor F^^, = \ e^yp^ Fp^. Using the definition of the 
dual field strength-tensor and Eqs. Q and (fTT|) yields for the magnetic four- 
current on the lattice 

jniA^) = ^^ef^^p^A^sinOp^lx) . (13) 

Evidently, this current is also conserved (A^jm,/^ = 0). The magnetic charge in 
the cube at position x is given by 

Qm{x) = a^jm,A{x) = ^ X] sin 6'p. (14) 

PG9cube(x) 

In this definition the monopole charge in each cube is not an integer multiple 
of g. Only for infinitely separated monopoles investigated in the classical limit, 
the charge inside a large surface enclosing a monopole is an integer multiple of g 
since sin^jj —>■ 6ij. This leads to the picture of a single cube yielding the main 
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contribution to the magnetic charge of the monopole, while the neighboring cubes 
contribute with a few percent. Since this charge distribution is independent of 
the lattice constant monopoles are still "point-like". This definition of monopole 
currents was investigated in detail in Ref. pi]. It was shown that the defini- 
tion (fT^ is necessary to get agreement with the dual superconductor picture of 
confinement and the dual London equation jT^. 

Using Eq. (fT^ we introduce a monopole's center of mass 

depending on the flux distribution around the monopole where r^^-* corresponds 
to space-like cubes with a positive contribution to the magnetic charge and cor- 
respondingly for f^~\ With definition (fT3j) it is possible to obtain positions of 
monopoles which are not restricted to integer values. 



5 Monopole Generation 

For the realization of a classical monopole-configuration on the lattice we discuss 
two possibilities: 

• Since a plaquette is part of the surface of four cubes, a non-vanishing pla- 
quette value corresponds to non-zero magnetic currents in these cubes. On 
the dual lattice the corresponding links assemble a plaquette. Following 
this line of thought, fixing the value of an a:?/-plaquette leads to a rotating 
current in the 2;t-plane. Further fixing of all x?/-plaquettes with given coor- 
dinates X, y, z and arbitrary t and closing them periodically at the boundary 
of the lattice yields two currents in t-direction but with opposite sign. For 
sufficiently large plaquette angles the two currents in t-direction represent 
a static monopole-antimonopole pair. Thus, arbitrary monopole configura- 
tions can be realized by fixing corresponding plaquettes appropriately. 

After specifying the constraints for a static monopole pair we minimize the 
action by local cooling. In fig. [T] the sum of the absolute values of the 
magnetic charges in one time slice 

WQmW = J2\Qm{r,t = to)\. (16) 
r 

in dependence on the angle of the plaquette separating the monopole pair 
is shown. 

Using definition |T0|) yields zero charge for plaquette angles in the inter- 
val (— vr, vr] followed by a discrete jump at ivr to 2g when a Dirac string 
penetrates the plaquette. Using definition (fT!^ this jump is spread out 
(see fig. [T]). This is caused by the fact that the charge of the monopole in 
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Figure 1: Monopole- charge in dependence on the plaquette-angle On- Besides 
the measured data and the result due to the definition of DeGrand- 
Toussaint an analytical approximation formula is plotted which is given 
by ||g™||/^7=^|^n-sin^n|. 



definition (fT3|l is distributed over several cubes, allowing for charge anni- 
hilation for a monopole pair at close distance. Furthermore, values for the 
approximation formula 

IIQ^II = -l^n-sin^nl (17) 

TT 

are plotted in fig. ^ This expression is obtained by considering two large 
closed surfaces each surrounding one monopole and assuming that all pla- 
quettes P, except the tuned one □, have small values 6p. 



• Constraints with local chemical potentials introduced at single cubes 
can be used to create monopoles at given positions. The potentials are 
introduced in the action as follows 

S, = 5^(1 - cos^p) + Y,^[C)J{C) , (18) 
p c 

with the plaquette angle Op, the chemical potential /^(C), the magnetic 
current J{C) and the sum C extending over all cubes. Fixing a chemical 
potential ±|/i| for a static monopole pair at a distance of one lattice spacing 
followed by local cooling preserves such a monopole-antimonopole pair for 
values > 2.3 only, see fig. 121 
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Figure 2: Absolute value of magnetic charge in dependence on the absolute value 
of the chemical potential introduced at neighboring cubes. 



One can observe that the charge for the monopole pair created at ~ 2.3 
is approximately 1.5(7 and increases continuously to a value near 2g. For 
smaller values of the chemical potential is not strong enough to keep 
apart a monopole pair at a distance of one against their magnetic attrac- 
tion. The monopoles approach each other with decreasing therefore the 
magnetic charge densities of monopoles and antimonopoles overlap leading 
to gradual annihilation. For increasing \^\ the monopoles are closer at the 
centers of the cubes and their overlap decreases. Since the charge of the 
monopole is distributed over more than a single cube a value of 2g cannot 
be attained for a monopole pair at distance one. In a strict sense, only 
infinitely separated monopoles obtain the value 2g. Furthermore, as soon 
as a monopole pair is created one would expect local minima of the action 
S^. Indeed this is true for |/i| > 2.3 as can be seen in fig. El 

In fig. 01 the action in dependence on the value On of the plaquette, which 
separates the monopole pair is shown for various values of One can clearly 
recognize local minima of for > 2.3. 

6 The Monopole- Antimonopole Potential in the Classical Limit 

The potential V oi a, static monopole pair in the classical limit corresponds to 
the action S ccV ~ T , since the kinetic energy T vanishes. Using the monopole 
definition (fTUI) due to DeGrand-Toussaint[Tn] the monopole charge is contained 
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Figure 3: S*^ according to Eq. ()18|1 for a monopole pair in neighboring cubes in 
dependence on the phase On of the plaquette separating the monopoles 
for various values of 



entirely in single cubes and the potential can be determined for integer distances 
only. Due to the electro-magnetic duality a monopole pair feels a 1/r-potential 
at large separations. 

We investigate the monopole-antimonopole potential not only for integer dis- 
tances but for arbitrary real ones. For this purpose monopoles must be movable 
inside cubes. This can be achieved by tuning of appropriate plaquettes. Fig. IH 
shows the result of such tuning. 

The chemical potential is set to /i = —2.5 at the cube centered at (—0.5, 0.5, 0.5), 
while the angle On of the plaquette with its center at (2, 0.5, 0.5) is tuned. Af- 
ter choosing appropriate values of Ou cooling generates an antimonopole near 
(—0.5, 0.5, 0.5) the center of the cube with non-vanishing /i. The position of the 
monopole depends strongly on the value of On- For 6'n > tt the plaquette is pierced 
by a Dirac-string and the position of the monopole is larger than two. With de- 
creasing On the monopole is continuously shifted towards the antimonopole. One 
can observe that for On = the monopole is shifted slightly to an x-value below 
2, which is due to magnetic attraction. Calculations for large separations with 
and without chemical potentials show that the Wilson part of the action does 
not differ. 

Using the above procedure it is possible to generate monopole pair configu- 
rations at arbitrary real distances. The potential of a static monopole pair in 
dependence on their distance yields fig. El 

The coordinate Xn of the tuned plaquette On varies from 1 to 5 in integer steps. 
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Figure 4: Relation between 9u and the positions of monopoles. The vertical 
straight line at x = 2 indicates the position of the tuned plaquette 
On- Introducing a chemical potential yU = —2.5 at a; = —0.5 and tun- 
ing of On followed by cooling leads to an antimonopole in the cube at 
X = —0.5 and a monopole in the vicinity of the tuned plaquette. 




Figure 5: The monopole-potential for arbitrary separations is shown. The dif- 
ferent symbols indicate the position Xn of the tuned plaquette which 
allows to vary the monopole position. 
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The values of Xn used to determine the potential are indicated in the figure. 
It is interesting to observe that the well-known 1/r-potential is superimposed 
by oscillations with a wave length of lattice spacing a. These oscillations are 
obviously due to the granularity of the lattice. Local minima are located at 
integer distances, while the maxima can be found at half-integer separations. 

This behavior of the potential at non-integer distances may lead to the well- 
known observation that some monopole loops survive the transition from the 
confined to the Coulomb phase. Refs. 7j and report strong indications that 
this phase transition is of first order. A vanishing of the two-state signal was de- 
tected in Refs. on a lattice with a topology homo-topic to an S^-sphere. 
This was explained by the fact that for such a topology monopole loops may 
contract and disappear without restriction. A discretization of an where the 
curvature is distributed almost homogeneously over the lattice favors a contrac- 
tion of arbitrary monopole loops and decreases the infiuence of meta-stabilities 
like those appearing in Fig. 

7 Conclusion and Outlook 

Magnetic monopoles defined by their charge distributions can take arbitrary po- 
sitions on a discrete space-time lattice. With this definition the potential can be 
calculated for arbitrary non-integer distances between monopoles. Local chemi- 
cal potentials provide an efficient tool to create monopoles on the lattice and to 
realize arbitrary monopole configurations. Magnetic monopoles feel a Coulomb 
potential in analogy to electric point-like charges, if they are located in the cen- 
ters of the lattice cubes. The above results demonstrate deviations from the 
1/r-behavior when the monopoles are shifted away from the centers. We find 
it interesting to investigate whether the meta-stabilities appearing in the poten- 
tial due to the resulting hysteresis are strong enough to infiuence the order of 
the phase transition. Methods to remove the infiuence of the granularity of the 
lattice will be a subject of our further investigations. 
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